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Abstract 

There has been great interest in applying the results of statistical me- 
chanics to single molecule experiements. Recent work has highlighted 
so-called non-equilibrium work-energy relations and Fluctuation Theo- 
rems which take on an equilibrium-like (time independent) form. Here 
I give a very simple heuristic example where an equilibrium result (the 
barometric law for colloidal particles) arises from theory describing the 
thermodynamically non-equilibrium phenomenon of a single colloidal par- 
ticle falling through solution due to gravity. This simple result arises 
from the fact that the particle, even while falling, is in mechanical equi- 
librium (gravitational force equal the viscous drag force) at every instant. 
The results are generalized by appeal to the central limit theorem. The 
resulting time independent equations that hold for thermodynamically 
non-equilibrium (and even non-stationary) processes offer great possibil- 
ities for rapid determination of thermodynamic parameters from single 
molecule experiments. 
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The equilibrium distribution for colloidal particles in dilute aqueous suspension 
follows a familiar barometric or exponential law ^ |2] 

c(/,) = C (0)ex P (-^) (1) 
k B T 

where c(h) is the concentration of particles at height h. For spherical particles of 
radius r the effective mass is m = 47rr 3 (p p — p w )/3 where p p and p w are the mass 
densities of the particle and of water, respectively. Eq. (1) is an equilibrium 
result that implicitly involves many particles so that the concentrations (particle 
densities) are well defined. 
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We can look at the situation of colloids from the very different perspective 
of a single particle falling through solution. The forces acting are gravity, mj, 
viscous drag jv ~ rjrv (where 77 w 10 -3 kg/(m s) is the viscosity of water), and 
a random thermal noise force the origin of which is the molecular movement 
of the water molecules. If we wait sufficiently long (about 771/7= 10~ 6 s for 
a micron sized particle) the particle will reach terminal velocity Wtcrm = mgH 
where the force of gravity is balanced by the viscous drag and there is no further 
acceleration. Then, taking as the origin of the coordinate system the center of 
mass of the particle at some time t — after terminal velocity is attained, the 
subsequent probability density function P(h, t) for the particle position can be 
written 

which is a solution of Fick's equation for diffusion with drift [3] 

dP(h,t) _ d 2 P{h,t) mg dP(h,t) 

dt ~ dh 2 + 7 dh ( ' 

where D is the diffusion coefficient. Eq. 2 describes a Gaussian distribution 
with mean /j, — (mg/j)t and variance a 2 = (h(t) 2 ) — (h(t)) 2 = 2Dt. Although 
it is more likely of course that the particle moves downward, at short times, 
there is a reasonable chance that thermal noise will cause the particle to be 
found slightly higher than where it started at t = 0. Since the distance from 
the center of the Gaussian to the original h = is [(mg/rft], the probability 

for a particle to be above its starting point at time t is ierfc [ ^ m fffi f ^ ] where 
erfc (x) is the complement of the error function. 

These "upward" trajectories have been termed "violations of the second law 
of thermodynamics on short time and small length scales" |1] because the en- 
tropy change is negative for trajectories where colloidal particles move up. The 
upward trajectories certainly do NOT demonstrate a violation of the more rele- 
vant Thomson (Lord Kelvin) formulation of the second law which states that no 
repeatable (cyclic) process can do work on the environment with the sole change 
being a decrease in the temperature of the system. The upward trajectories can 
however be exploited in a non-isotropic system to allow the random input of 
energy to drive directed motion by a Brownian motor mechanism 

We can calculate exactly the ratio of the probability density for the particle 
to "fall up" to +h to the probability density to fall down to —h from Eq. 2 to 
be 

P(h,t) _ exp( Z^) (4 ) 



P(-h,t) ^ 1 D 

Remarkably, time has disappeared altogether in this ratio, and by inserting 
Einstein's relation |B] 7D = ksT we regain the "equilibrium" barometric law 
Eg. l! An analogous relation for motion of an overdamped particle in an 
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arbitrary potential was derived in a more general context by Bier et al. \7\ 
using Onsager's [SJ thermodynamic action approach. 
Since the probability density function is normalized 



/-(-oo r+oo 
P(h)dh = I P{-h)dh = 1 

-oo J — oo 



we also have 



(«pS)>=/ 00 e XP Ap(h)dh=l (5) 
k BT k B T 

The quantity Wdiss = —rngh is the work dissipated when a particle falls a 
distance h. 

Many stochastic processes obey a Gaussian distribution due to the central 
limit theorem 9 which states that the probability density function for any 
variable that can be expressed as the sum of many small random quantities ap- 
proaches a Gaussian distribution irrespective of the microscopic laws governing 
the individual small quantities so long as the mean and variance of the small 
quantities are finite. Consequently we expect equilibrium-like (time indepen- 
dent) relations to hold for a wide range of processes. 

Consider an arbitrary macroscopic variable q, often termed a generalized 
displacement, and the time derivative of its average q which is often called 
a generalized velocity or flux. If q fluctuates due to molecular noise in the 
environment the distribution for q will be given by 

P(9,*)= ,™ 1 (6) 
ViirDt 

where the variance a = (q(t) 2 ) - (q(t)) 2 = 2Dt and we take q(0) = 0. As before 
P(q, t) = P(-q, t) exp [qt/D] and 

(exp [<#/£]> = !• (7) 

Further, if the generalized velocity q is related linearly to a generalized force X 
such that q — X/j (as is the case after terminal velocity is reached) we have 
P{q,t) = P(-q,t) exp [qX/faD)]. Finally, with identification of 7D = k B T (a 
fluctuation-dissipation relation) we have 

/>(,/i -exp[(WW(fc B T)] (8) 



P(-q) 



where as before the dissipated work Wdiss is the product of the generalized force 
and displacement —qX. Using the thermodynamic relation that the dissipated 
work is the difference between the total work and the stored free energy Wdiss = 
W — AG we have the equality 



(exp (W/(k B T))} = exp(AG/k B T) 



(9) 
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Relations analogous to equations (8) and (9) (termed non-equilibrium fluctuation- 
dissipation theorems (FDT)) derived from a different perspective were published 
twenty five years ago by Bochkov and Kozovlev |10| . 

Interest in the non-equilibrium FDT has recently been rekindled by the 
theoretical work of Jarzynski [11] . Bier et al. [Jj, Crooks ^2]> and others. 
This work has made important progress by relating experimental observables 
from non-equilibrium experiments to thermodynamic parameters such as free 
energy. Hummer and Szabo pointed out that the relations are implicit 
in the Feynman-Kac path integral theorem, and demonstrated how the non- 
equilibrium work energy relations and fluctuation theorems can be used to in- 
terpret single molecule experiments. The relations (8) and (9) have been exper- 
imentally tested by Bustamante and colleagues |141 1151 HB] who showed that a 
relatively small number of single molecule experiments can be used to obtain 
very good estimates of the Equilibrium Free Energy profiles for the system. 

It is often stated that many single molecule experiments are carried out far 
from equilibrium. Indeed, the 10-100 nN force typical for an AFM pulling ex- 
periment is enormous compared to the, say, 10-20 pN force associated with even 
strong molecular motors [T7j and the rate of change of the force - up to several 
hundred nN per second - seems very large. However, when we compare the rate 
of change of the force (dF/dt ~ 100 nN/s) with the ratio of the characteris- 
tic force JB] -Fchar = V 2 I P ~ 10 _9 N to the characteristic relaxation time for a 
nanometer object in solution i c h ar = P r 2 /rj oc m/7 ~ 10~ 12 s (-F c har/tchar ~ 10 3 
N/s) we see that the system is in fact very close to mechanical equilibrium at 
every instant in time. The characteristic quantities are written in terms of the 
particle size r ~ 10~ 9 m, density p ~ 10 3 fc<?/m 3 , and the viscosity of the so- 
lution 77 ~ 10~ 3 fc<?/(ms). Nevertheless, the system is far from thermodynamic 
equilibrium as it continues to move, and to dissipate energy in the environment. 

In nanoscale physics, it is essential to remember that "equilibrium" does 
not have a single unambiguous meaning. The thermodynamic non- equilibrium 
aspect of a typical experimental system shows up as a mean energy dissipa- 
tion rate that is not zero. The mechanical equilibrium aspect shows up in the 
Gaussian distribution of fluctuations about the mean dissipated work (or mean 
extension) over many realizations. The MECHANICAL equilibrium allows a 
simple description of the distribution of the fluctuations of quantities about 
their mean values in THERMODYNAMICALLY far from equilibrium situa- 
tions. The closeness to mechanical equilibrium is ultimately the explanation 
for the unreasonable effectiveness of equilibrium theory for interpreting experi- 
ments which, from our perspective (but not from that of the single molecule!) 
appear to be "far from equilibrium" . 

Time independent equations that hold for non-equilibrium (and even non- 
stationary) processes offer great possibilities for rapid determination of thermo- 
dynamic parameters from single molecule experiments. The present paper is 
intended to provide a simple background for understanding these equilibrium- 
like relations, and shine light on their historical origins to facilitate further 
investigation. 
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